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Dynamical structure factors of the S = 1 bond-alternating spin chains in the dimer phase are 
calculated at finite temperature, using the pair dynamical correlated-effective-field approximation. 
At T = 0, the delta-function-type peak of the one-magnon mode appears. When temperature is 
increased, such a sharp peak is broadened and the additional weak peak caused by the excitation 
from the triplet state to the quintet state emerges in the higher energy region. The results are 
£^ . discussed in comparison with those obtained by the exact diagonalization method. 
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I. INTRODUCTION 

, There has been a continued interest in S — 1 bond-alternating Heisenberg chains both theoretically and experi- 
mentally. It was shown that a quantum phase transition between the Haldane phase and the dimer phase takes place 
depending on the bond-alternating ratio P, Q . The quantum phase transition of the system was studied by including 
i the effects of a single-ion anisotropy. The phase diagram for the parameters of the bond-alternating ratio (a) and the 
single-ion anisotropy (D) was obtained 0,0,0, where the bond-alternating ratio is 1 : a. Along the phase boundary 
of the Haldane-gap and dimer phases, the system becomes gapless. 

Recent progress on material synthesis has made it possible to study the elementary excitation as well as thermo- 
dynamic properties of S = 1 bond-alternating Heisenberg chains experimentally. Furthermore, dynamical properties 
have been investigated by inelastic neutron-scattering experiments for a typical anisotropic dimer-phase compound 
Ni(C 9 D24N4)(N0 2 )C104 (abbreviated to NTENP) [gj. It was shown that spin dynamics of NTENP in transverse 
magnetic fields show noticeable different features as compared with those of an anisotropic Haldane-gap compound 
| Ni(C5Di4N2)2N3(PFg) Q- The experimental findings have been analyzed on the basis of the numerical diagonalization 
. calculation and explained in viewpoint of the field dependence of the excitation continuum [j| . 
in ' Synthesized compounds for S = 1 bond-alternating Heisenberg chains so far are in the dimer phase. NTENP is 
located relatively close to the gapless line. Some compounds are located far away from the gapless line. For example, 
the parameters of [Ni2(dpt)2(^ — ox)(^ — N3)](PF@)[dpt = bis(3 — aminopropyl) amine, ox = C2O2] (abbreviated to 
NDOAP) are evaluated as a — 0.1 and D = and the system is situated close to the isolated dimer system 0. 
Dynamical structure factors (DSF) of S = 1 bond-alternating Heisenber g ch ains at zero temperature were calculated, 
using the continued fraction method based on the Lanczos algorithm [8lllf)j|. Field dependence of spectral intensities 
I"-; . of the magnon isolated mode and the excitation continuum was investigated in the Haldane-gap and the dimer phases. 
However, temperature dependence of the dynamical structure factor has not yet been fully investigated. It may be 
desirable to investigate dynamical properties of S = 1 bond-alternating Heisenberg chains at finite temperatures. 

In this Letter, we calculate dynamical structure factors at finite temperatures in the dimer phase, using a pair 
dynamical correlated-effective-field approximation (Pair-DCEFA) Pair-DCEFA was successfully applied to cal- 
culate the dynamical structure factor as well as the static susceptibility of the S = 1/2 bond- alternating Heisenberg 
chain at finite temperatures |ll|. In Section 2, we derive the effective Hamiltonian and formulate the expressions for 
the static susceptibility and the DSF, by using Pair-DCEFA. In Section 3, we first show the numerical results for 
the DSF at zero temperature and the temperature dependence of the susceptibility. The results are compared with 
those obtained by numerical diagonalization method. We then show the results for the DSF at finite temperatures. 
$J] ■ Characteristics of the isolated modes are discussed. Brief summary is given in the last section. 
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II. MODEL AND METHOD 

Let us consider the S = 1 bond-alternating Heisenberg chain described by the following Hamiltonian, 

H = jj(SySy+aSi.y^) (J > 0), (1) 



where Si. 1 (Si, 2) denotes the operator of the spin on the left side (right side) in the unit cell i. 
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In Pair-DCEFA, we first consider the singlet, triplet, and quintet states of an isolated spin pair on the strong bond. 
To take into account the effects of the interaction between spin pairs, we decouple the second term of the Hamiltonian 
JU and then obtain the effective Hamiltonian in Pair-DCEFA as 



H cS = J^2 S i,i ■ S *,2 + ajY, [Si,i • «Si_i, a > - A(S M )) + S i;2 • «S i+M ) - A<S,, a » 



(2) 



where (Sj U ) (y = 1, 2) denotes the spontaneous moment and A is a correlation parameter determined by the self- 
consistency equation derived from the fluctuation-dissipation theorem. We confine ourselves to consider the param- 
agnetic state. On the basis of the effective Hamiltonian, the self-consistency equation for A is obtained as 



coth 
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coth 
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diTiL)+(q) 
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d 2 J 2 
hu- (q) 
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(3) 



where N is the total number of unit cells, j3 = l/k B T, d x = (4/3)(2Ap + 5Apx), and d 2 = (4/3)(8Ap + 5Api) 
with Ap = Ps — Pt and Api = p t — p q . Here p s (= exp(2(3J)/Z), p t (= exp(/3J)/Z), and p q {= exp(—(3J)/Z) are the 
contributions from the singlet, triplet, and quintet states with Z = exp(2/3J) + 3exp(/3J) + 5exp(— (3 J). In thermal 
equilibrium, the excitation from the singlet state to the triplet state and the excitation from the triplet state to the 
quintet state take place. When T — > 0, we find that p a — > 1, pt — > 0, and p q — > 0. 
The dispersion relation of the one-magnon mode is given by 



J 



LU±(q) = —j= (5 - ad\{\ + cos[g(c + c')]}) 
V 2 L 



± J (5-adi{A + cos[ 9 (c + c')]}) -4(4 - ad 2 { A + cos[ g (c + c')]}) 



(4) 



where w+(q) and 0J-(q) denote the dispersion relations of the excitations from the singlet state to the triplet state 
and from the triplet state to the quintet state, and c and c' are the lengths between the neighboring spins along the 
strong bond and the weak bond, respectively. 

By using A and to±(q) thus obtained, the static susceptibility x(T) and the dynamical structure factor S(q, w) in 
Pair-DCEFA are expressed as 



X(T) = 



2(g^B) 2 (pt + 5 Pq ) 

1 + (A - l)aJ(p t + 5p q ) ' 



(5) 



o, , 2J 2[l-coB(gc)] r d^Tiu+iq)} 2 - d 2 J 2 

o{q,uj) 



1 _ e -/3Ku, [hui+(q)] 2 - [Tiw-(q)] 2 I ncJ + (q) [u - u+{q)] 2 + 5 2 

d^fiu^q)] 2 - d 2 J 2 S 



-}■ (6) 



hu_{q) [Lj-Lj-(q)] 2 +S 2 . 

In the following, we set c = c' = 1, fee = h = 1, and S = 5.0 x 1CP 2 . The energy is measured in units of J. 



III. RESULTS 



We solve Eqs. (jSJ and (@J numerically to obtain A and o-i±(q) in given T for fixed a. The results for A are shown 
in Fig. ^ In |a| > 0.6, the convergence of the self-consistent procedure becomes worse at finite temperatures. 
Such ineffectiveness of the self-consistent treatment may relate to the quantum phase transition around a ~ 0.6 
[l2l IT3L IT4I IT5I IT?! lri| . The behavior of A in this system is in contrast with that in the S = 1/2 bond-alternating spin 
chain, where temperature dependence of A is obtained even at |A| = 1 [TT| . The different feature of A may represent 
characteristics of spin correlations in both systems. 

We next show the dispersion curve o-i-(q) at T = 0. In Fig. [21 the results are expressed in the extended zone 
scheme. The dispersion curves obtained by Pair-DCEFA are expressed by the solid lines for a = —0.3, 0.1, 0.3, and 
0.4. As a reference, the excitation energies obtained by the exact diagonalization method for 20 spin systems [nj are 
plotted by dots. For a = 0.1 and 0.3 the agreements between both results are excellent, while for a = —0.3 and 0.4 
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FIG. 2: The dispersion curves uj-(q) at T = for a = -0.3,0.1,0.3, and 0.4. The results obtained by Pair-DCEFA are 
expressed by the solid lines and those obtained by the exact diagonalization method for 20 spin systems are plotted by dots. 



the deviations become noticeable in q ~ and it. Thus, Pair-DCEFA may be effective in —0.3 < a < 0.3 for the 
5 = 1 bond-alternating Heisenberg chain. 

Temperature dependence of the static susceptibility x(^) i s shown in Fig. The results for a = —0.1,0.1, and 
0.3 are compared with those obtained by the exact diagonalization method for 8 spin systems. Note that %(T) by 
the exact diagonalization method exhibits little size dependence around 8 spin systems. The agreement between both 
results becomes better in low temperatures. 

Using numerical results for A and u>±(q), we calculate dynamical structure factors at finite temperatures. The 
results are shown in Fig. 0] According to the finite-size analysis for the S(q, ui) obtained by the exact diagonalization 
method, the lowest excited states in a — —0.1, 0.1, and 0.3 form the one-magnon mode in < q < n [Icj . 

Since we turn our attention to the behavior of the isolated mode located around u ~ 1.0 and 2.0, we set T < 1.0. 
At T — the intensity in a = —0.1 takes the maximum at q — 0.957T, while in a — 0.1 and 0.3 the intensity at T = 
takes the maximum at q = n. We have confirmed that the wave number of the maximum intensity at T = shifts 
towards 7r/2 , as a decreases in a < 0. The results are consistent with those obtained by the exact diagonalization 
method [Tg. 

When temperature increases, the intensity of the one-magnon mode around to ~ 1.0, which corresponds to the 
singlet-triplet excitation, is reduced. At finite temperatures, the additional mode with weak intensity emerges around 
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FIG. 3: The static susceptibility x(T) for a = —0.1,0.1, and 0.3. The results obtained by Pair-DCEFA are expressed by the 
solid lines and those obtained by the exact diagonalization (ED) method are expressed by the broken lines. 



u> ~ 2.0. At T = 0.8, its intensity takes the maximum at q ~ tt in a = 0.1 and at q = 0.787T in a = 0.3, respectively, 
whereas the intensity of the singlet-triplet mode in a > takes the maximum at q = tt even at finite temperatures. 
The excitation mode around u) ~ 2.0 is caused by the triplet-quintet excitation in thermal equilibrium, which vanishes 
at T = 0, because pt — > and p q — > as T — > 0. 

As shown in Fig. |3 the excitation continuum is located around u> ~ 2.0 in a = 0.3 [ljj. Accordingly, the triplet- 
quintet mode becomes unstable in the continuum and disappears. In a = 0.1, by contrast, the excitation continuum 
is located in the higher energy region |10|. Therefore, the triplet-quintet mode may be stable. Recently, it was 
determined experimentally that NDOAP is well described by the isotropic S = 1 bond-alternating Heisenberg chain 
with a ~ 0.1 9]. When inelastic neutron-scattering experiments are performed on NDOAP, such a triplet-quintet 
mode with weak intensity may be observed around w ~ 2.0 in addition to the noticeable singlet-triplet mode around 
lo ~ 1.0. According to our results, the ratio of the intensity of the higher-energy mode to that of the lower-energy 
mode at q = it is evaluated to be 22% in T = 0.8. 



IV. SUMMARY 



We have investigated dynamical properties of S = 1 bond-alternating spin chains in the dimer phase at finite 
temperatures, using Pair-DCEFA. At low temperatures, the delta-function-type peak of the excitation from the 
singlet to triplet states appears in w ~ 1. When temperature increases, the additional weak peak caused by the 
excitation from the triplet to quintet states emerges in oj ~ 2. This additional mode may be observed in the system 
with small a such as NDOAP. 
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